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Abstract 


Finite  deformation  in  the  crack-tip  zone  of  plastic  deformation  is 
investigated  for  Mode-I  opening  of  a  crack  in  a  thin  sheet  of  elasto-plastic 
material.  The  material  obeys  the  von  Mises  yield  criterion  in  the  true 
stresses,  and  the  stretching  tensor  satisfies  a  flow  law  of  the  Prandtl- 
Reuss  type.  I^ompressibility  and  a  state  of  generalized  plane  stress  are 
assumed.  It  is  assumed  that  linearized  elasticity  applies  outside  the  zone 
of  plastic  deformation.  On  the  crack  line  between  the  crack  tip  and  the 
elastic-plastic  boundary,  two  distinct  regions  have  been  recognized;  the 
near  tip  zone  and  the  intermediate  region.  In  the  near  tip  zone  the  fields 
are  controlled  by  the  radius  of  curvature  of  the  blunted  crack  tip.  Here 
the  stress  field  has  been  approximated  by  classical  plane  stress  results. 

It  has  been  assumed  that  the  crack-line  stresses  may  be  taken  as  uniform 
in  the  intermediate  region.  In  each  region,  deformation  variables  have 
been  determined  by  the  use  of  the  constitutive  relations,  and  the  results 
have  been  matched  to  the  corresponding  quantities  in  the  neighboring  region(s) 
In  thi;  manner  expressions  have  been  constructed  for  the  deformation  gradients 
on  the  crack  line,  in  terms  of  the  distance  to  the  crack  tip  in  the  deformed 
configuration,  the  yield  stress  in  shear  and  the  stress-intensity  factor  of 
linear  elastic  fracture  mechanics. 
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Introduction 

It  is  difficult  to  include  the  effects  of  deformation  induced  geometry 
changes  in  an  analysis  of  stresses  and  deformations  near  a  crack  tip.  For 
plane-strain  conditions.  Rice  and  Johnson  (1970)  considered  the  modification 
of  elasto-plastic  crack-tip  stress  and  strain  states  due  to  crack  blunting. 
McMeeking  (1977)  and  Rice  et  al  (1979)  used  finite  element  methods  to  include 
finite  deformation  in  their  investigations  of  crack-tip  opening  in  elastic- 
plastic  materials.  Finite  deformations  near  a  stationary  crack  tip  in  an 
elastic  solid  have  been  investigated  by  Knowles  (1977),  Knowles  and  Sternberg 
(1980,  1981),  and  Lo  (1977),  who  used  analytical  methods. 

In  this  paper  we  present  an  approximate  analysis  of  the  effect  of  finite 
deformation  on  the  crack-line  fields  for  Mode-I  opening  of  a  crack  in  a  thin 
sheet  of  elasto-plastic  material.  Away  from  the  crack  tips  the  strains  are 
assumed  to  be  infinitesimal,  and  the  material  behaves  according  to  the  theory 
of  linear  elasticity.  The  zones  of  plastic  deformation  near  the  crack  tips 
are  small  as  compared  to  the  crack  length  (small-scale  yielding).  Finite 
deformation  is,  however,  taken  into  account  in  the  near-tip  plastic 
deformation.  The  constitutive  model  for  the  elastic-plastic  deformations  is 
a  finite  deformation  version  of  the  elastic  perfectly-plastic  solid,  with  the 
Prandtl-Reuss  flow  rule  and  the  von  Mises  yield  condition.  It  is  assumed 
that  the  material  is  imcompressible ,  and  that  the  stress  state  in  the  sheet 
may  be  approximated  by  generalized  plane  stress. 

We  consider  a  crack  which,  in  the  reference  state,  is  defined  by  <  0, 

=  0.  The  crack  faces  remain  free  of  surface  tractions.  In  the  deformed 
configuration  the  stress-equations  of  equilibrium  have  the  same  form  as  for 
the  corresponding  small  strain  theory.  It  is  assumed  that  stressing  of  the 
sheet  gives  rise  to  a  blunted  crack  tip  profile  of  smoothly  varying  radius 
of  curvature,  whose  maximum  value  is  reached  on  the  crack  line.  We  focus  on 


the  fields  near  the  crack  line,  and  we  consider  a  near-tip  zone  and  an 
intermediate  region  in  between  the  near-tip  zone  and  the  elastic-plastic 
boundary.  In  the  near-tip  zone  the  fields  of  stress  and  deformation  should 
be  greatly  affected  by  the  radius  of  curvature  of  the  blunted  crack  tip.  In 
this  zone  we  use  a  classical  stress  field  that  has  been  applied  previously 
to  a  sheet  containing  a  circular  hole  and  for  a  sheet  with  an  edge  notch, 
see  Kachanov  (1971).  This  field  is  axially  symmetric.  In  the  intermediate  region 
we  assume  uniform  crack-line  stresses.  It  is  to  be  expected  that  the  actual 
stresses  near  the  crack  line  will  develop  a  dependence  on  the  polar  angle 
in  the  transition  from  the  essentially  axisymmetric  field  very  close  to  the 
blunted  tip  to  the  uniform  field  in  the  intermediate  region.  For  simplicity 
we  have , however ,  assumed  that  the  axisymmetric  field  remains  applicable  until 
it  can  be  matched  to  the  uniform  stresses  in  the  intermediate  region.  The 
choice  of  the  stress-field  in  the  intermediate  region  is  shown  to  be  consistent 
with  a  slip-line  field  proposed  by  Thomason  (1979). 

When  the  stress  fields  have  been  chosen,  the  crack  line  deformation  can 
be  analyzed  by  the  use  of  the  constitutive  relations.  The  deformation  gradients 
can  be  solved  rigorously  in  the  near-tip  zone  in  terms  of  a  single  function 
of  the  crack-tip  radius,  whose  form  is  subsequently  determined  by  matching 
deformation  variables  to  corresponding  ones  in  the  intermediate  region.  In 
the  latter  region  the  crack-line  deformation  is  analyzed  by  following  the 
method  discussed  by  Achenbach  and  Dunayevsky  (1984)  and  Achenbach  and  Li  (1984). 
The  solution  is  completed  by  matching  the  intermediate-region  results  to  the 
deformation  gradients  at  the  elastic-plastic  boundary. 
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The  principal  results  of  the  paper  are  expressions  for  the  deformation 
gradients  in  terms  of  the  distance  to  the  crack  tip  in  the  deformed  con¬ 
figuration,  the  yield  stress  in  shear  and  the  stress  intensity  factor  of 
linear  elastic  fracture  mechanics.  It  is  shown  that  is  singular 

at  the  crack  tip,  while  Sx^/SX^^  vanishes  at  that  point. 
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1 .  Governing  Equations 

The  reference  and  current  positions  are  defined  by  postion  vectors 
X  and  X,  respectively.  In  the  current  configuration  the  equilibrium  equation 


divT  =  0  , 


(1.1) 


where  T  is  the  Cauchy  stress.  In  this  paper,  the  spatial  derivatives  are 
with  respect  to  x  unless  stated  otherwise.  The  constitutive  model  that  is 
being  considered  is  a  large  deformation  version  of  the  elastic  perfectly- 
plastic  solid,  with  the  Prandtl-Reuss  flow  rule  and  the  von  Mises  yield 
criterion.  In  addition  it  is  assumed  that  the  material  is  incompressible. 

The  yield  condition  is  represented  by 


T' -T'  =  2k^  , 


(1.2) 


where  k  is  the  yield  stress  in  pure  shear  and  T'  is  the  stress  deviator 
defined  by 


T'  =  T  -  ^  I  trT  , 


(1.3) 


Here  I  and  tr  stand  for  the  unit  tensor  and  the  trace,  respectively.  When 
(1.2)  is  satisfied,  we  use  an  incremental  constitutive  equation  of  the  type 


D  =  ^  (T  +  pi)  +  AT' 

-  /'ll  -  ‘  ~  - 


(1.4) 


The  elastic  part  of  D  is  represented  by 


(T  +  pi)  . 


(1.5) 
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Here  D  represents  the  stretching  tensor,  which  is  defined  by 

D  =  i  [Vv  +  (Vv)”^]  ,  (1.6) 

where  v  is  the  particle  velocity: 

3x 

V  =  X  =  —  (X(x,t),t)  .  (1.7) 

The  dot  symbol  denotes  the  material  time  derivative.  In  (1.4)  and  (1.5) 
p  defines  the  indeterminate  pressure,  A  a  positive  multiplier,  and  the  ° 
symbol  defines  the  Jaumann  rate: 

O  • 

T  =  T  -  WT  +  TW  ,  (1.8) 

where 

W  =  y  [Vv  -  (Vv)'^]  .  (1.9) 

It  should  be  noted  that  there  is  some  controversy  over  the  choice  of 
stress  rate  in  rate-type  constitutive  equations,  see  e.g.  Lee  et  al  (1983)  and 
Dienes  (1979).  If  one  wishes  to  rewrite  a  small  deformation  constitutive 
equation  of  rate  type  into  its  large  deformation  counterpart  just  by  replacing 
the  small  deformation  stress  rate  by  an  objective  stress  rate,  it  is, to  the  present 
authors'  opinion,  generally  inadvisable  to  use  the  Jaumann  rate.  This  is 
because  W  is  not  a  measure  of  the  rate  of  rotation  of  a  material  point.  Indeed, 
it  is  not  difficult  to  find  a  deformation  gradient  F  which  is  symmetric  positive 
definite  (rotation  -  free  with  regard  to  polar  decomposition) ,  but  with  non-zero 
W.  Actually 

F  =  Q(t)rQ^(t) 


(1.10) 


provides  such  an  example,  where  Q(t)  Is  a  time  dependent  orthogonal  tensor, 
aid  r(7^  yI>  Y  =  const,)  is  a  symmetric  positive  definite  constant  tensor. 

From  the  physical  point  of  view,  it  seems  that  in  the  constitutive  equations 
(1.4)  and  (1.5)  ,  Dienes'  rate.  Dienes  (1979),  is  a  more  reasonable  choice 
than  Jau,mann's  rate.  In  the  present  analysis,  however,  the  use  of  Jaumann's 
rate  is  acceptable  because  the  principal  axes  of  stress  do  not  rotate 
significantly.  In  addition,  Jaumann's  rate  is  much  easier  to  handle  than 
other  rates.  These  are  the  reasons  for  using  (1.4)  and  (1.5)  here. 

The  pressure  p  can  be  eliminated  from  (1.4)  and  (1.5)  by  using  the 
condition  of  incompressibility 


div  V  =  0  .  (1 

By  taking  the  trace  of  (1.4)  or  (1.5)  we  obtain 

p  =  -  i  tr  T  ,  (1 

which  in  turn  reduces  (1.4)  to 

D  =  ^  (T  -  j  I  tr  T)  +  AT'  ,  (1 

In  the  present  paper  we  will  discuss  the  deformations  of  a  thin  cracked 
sheet  of  uniform  initial  thickness  h.  As  in  the  corresponding  small 
deformation  theory,  one  can  facilitate  the  analysis  by  using  the  generalized 
plane-stress  assumption.  To  see  the  implication  of  this  assumption  we 
introduce  a  Cartesian  coordinate  system  with  its  =  0  plane  coinciding 
with  the  mid-plane  of  the  undeformed  sheet.  Under  the  assumptions  that  the 
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deformation  is  symmetric  with  respect  to  the  mid-plane,  that  the  1st  Piola 
Kirchhoff  stress  does  not  vary  significantly  throughout  the  thickness,  and 
that  the  planer  faces  are  free  of  loads,  Knowles  and  Sternberg  (1983)  showed 
that 


T  „  =  T-  =  0  , 
a3  3a 


F-  =  F  -  =  0  . 
3a  a3 


T  =0 
33 


(1 


hold  on  =  0,  where  Greek  indices  stand  for  either  1  or  2 .  If  (1.14 a) 
holds  not  just  on  =  0  but  in  its  vicinity,  we  may  reduce  (1.1)  to  its 


2D  version,  i . e. , 


This  equation  and  (1.2)  constitute  the  governing  equations  for  the  stress 
components.  It  is  noted  that  the  same  equations  govern  the  stress  dis¬ 
tributions  in  the  theory  of  small-strain  plane-stress  plasticity.  Therefore, 
we  may  utilize  some  known  results  from  that  theory.  Another  implication  of 
the  plane  stress  assumption  is  that 


V  =  V  „  =  0, 

3, a  a, 3 


(x^  =  0). 


With  this  result  and  the  constitutive  equations  (1.13),  we  can  determine 
the  velocity  field  in  the  mid-plane. 


( 


.  l4a,b,c) 


.15) 


2.  Crack-Line  Stresses 


We  start  with  the  observation  that  the  equilibrium  equation  (1.15)  and 
the  yield  condition  (1.2)  have  the  same  form  as  for  small  deformation  theory, 
except  that  in  the  present  context  the  stress  is  referred  to  the  x  position 
(after  deformation) .  Therefore  the  classical  plane  stress  analysis  for 
stresses  applies  to  the  present  problem,  but  with  the  boundary  conditions 
imposed  on  deformed  boundaries. 

It  is  assumed  that  the  crack  tip  will  blunt  into  a  curve  S  which  has  a 
root  radius  whose  magnitude  attains  its  maximum  value  a  on  the  crack  line. 

For  small  polar  angles,  we  may  then  approximate  S  by  a  segment  of  a  circle 
C  of  radius  a. 

At  this  point,  we  recall  some  formulas  from  the  classical  theory  of  plane 
stress  plasticity,  see  e.g.  Kachanov  (1971),  pp.  262-272.  A  set  of  principal 
stresses 

T,  =  2kcos(a)--^),  T.  =  2kcos(w+-^),  (0<w<it),  (2.1a,b) 

J.  b  I  D  —  — 

satisfies  the  yield  condition  (1.2),  where  u)  is  a  parameter.  When 
i/b  <  ‘ji  <  the  equilibrium  equation  becomes  hyperbolic.  The  equations 

for  the  characteristics  and  the  corresponding  integrals  are: 
a  -  characteristics 

dx2 

-  =  tan(d)-v)  ,  2  -  >1>  =  const  ,  (2.2a,b) 

d^l 

D  -  characteristics 
dx., 

=  tan(0+v)  , 

dx. 


+  ^)  =  const  , 


(2.3a,b) 


9 


where  (}>  is  the  angle  between  the  and  directions, 


and 


-  Y  arccosC-"^^*^)  , 


(2.4) 


tu 

--U. 


(3-4cos^C)*'^ 

sinC 


dC. 


(2.5) 


With  regard  to  the  crack-line  stresses,  one  would  expect  that  the 
stress  right  near  the  tip  is  affected  by  the  traction  free  boundary  S, 
while  the  stress  away  from  the  tip  (but  still  within  the  plastic  zone) 
will  not  be  all  that  different  in  form  from  those  for  the  corresponding 
small-strain  crack  problems.  This  would  suggest  that  we  should  approximate 
the  crack-line  stresses  by  using  two  fields:  one  for  the  immediate  vicinity 
of  the  crack  tip,  and  another  for  the  intermediate  region  between  the  near¬ 
tip  zone  and  the  elastic-plastic  boundary. 

2 . 1  Stresses  in  the  near-tip  zone 

Let  A  and  B  be  points  on  S  between  which  C  and  S  do  not  differ  significantly, 
see  Fig.  1.  Also,  let  D  be  a  domain  surrounded  by  the  arc  AB,  and  two 
characteristics  defined  by  (2.3)  and  (2.2),  emanating  from  A  and  B.  It  is 
reasonable  to  approximate  S  between  A  and  B  by  C  .  Then  the  stresses  in  D 
are  found  as  solutions  of  the  equilibrium  equation  (1.15)  and  the  yield 
condition  (1.2)  subject  to 


T  =  0  , 
r 


=  0 


on  r  =  a 


(2.6) 


where  T^,  T,,  and  T^ stand  for  the  components  of  T  referred  to  the  polar 


I  T>  ■  'i  V  ■  .1  ■  V  fiVi  tu... 


coordinate  system  (r,0)  which  has  its  origin  at  the  center  of  C.  Due  to  the 
hyperbolicity  of  the  governing  equation,  vanishes  identically  within  V,  and 
and  Tg  show  axisyinmetry .  The  solution  to  this  problem  appears  in  standard 
textbooks  of  plasticity  in  connection  with  the  extension  of  a  sheet  with  a 
circular  hole.  We  record  the  following  results  from  Kachanov  (1971); 


/r  \2  ^  /3  ^v'3(TT/3-aj) 

\3L  J  2sin(jj 


{2.7a) 


T  =  2k  cos  (u  +  — )  ,  T  =  2k  cos  (w  -  ^)  (0  <  uj  <  ^/3) 

r  0  0  0 


(2.7b,c) 


2 . 2  Crack-line  stresses  in  the  Intermediate  region 

Two  Mode-I  near-tip  stress  distributions  have  been  proposed  for  the 
small-deformation  plane-stress  case  in  an  elastic  perf ectly-plastic  material. 
Hutchinson  (1968)  proposed  a  stress  field  which  corresponds  to  the  system  of 
characteristic  curves  shown  in  Fig.  2.  He  also  showed  that  this  distribution 
can  be  obtained  as  the  limit  for  a  certain  non-linear  elastic  material  as 
the  stress-strain  curve  becomes  flat.  In  the  loading  zone  ahead  of  the  crack 
tip  the  stresses  are 

T^j^  =  k  cos^G,  T^^  =  k(2cos^0  +  3sin"acos0),  T^^  =  -k  sin^G.  (2.8a,b,c) 

The  crack  line  stresses  are  T^^  =  k,  T^,,  =  2k  and  Tj^^  =  0.  The  latter 
results  also  follow  from  (2.1),  since  it  can  be  deduced  from  (2.2)-(2.4)  and  Fig.  2 
tliat  the  crack  line  is  a  characteristic  curve  with  =  tj/6. 

Thomason  (1979)  proposed  a  class  of  stress  fields  which  include  the 
Hutchinson  field  as  a  special  case.  The  system  of  characteristic  curves  is 


(5.11) 


F(a)  =  —  a  . 

Substitution  of  (5.11)  into  (3.28)  -  (3.30)  completes  the  desired 

relations  on  the  crack  line  in  the  near-tip  zone,  a  <  <  Ka.  In  the 

intermediate  region,  xa  <  x.  <  x  ,  we  find  from  (5.11),  (5.7)  and  (3.28) 

1  P 


X  =  (x  -  v^a/x) 
1  1 


(5.12) 


By  using  (4.13),  (5.2),  (5.1)  and  the  previously  computed  result  for 


C(t),  the  expression  for  v  becomes 

Z  9  Z 

*  _  ^''3a  1+6  //3  ^  \  2/3a  6 

^2,2  3  2  1-6  \K  ^  ^o)  2  l-S 


(5.13) 


For  <a  <  x^  <  X  ,  (4.14)  then  takes  the  form 
1  P 

Sx-  *  -qx  [l/a-l/G*(x  )] 

-  F(x  )ta/G*(K  jl”  e  ° 
d.X2  o  o 


(5.14) 


where 


3(1+6)  -  2/3<^6 
<“*(1-6) 


(5.15) 


v'3(l+6) 
<'  (1-6) 


(5.16) 


* 

G  (x  ) 
o 


G(x  ) 

/  °  . 


(5.17) 


The  function  G  (x  )  mav  be  selected  in  a  convenient  manner  since  C'+y,  )  is  arbitrary, 

O  -  O' 

see  also  the  footnote  under  Eq.(4.12). 
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X.  =  F(x  /-)  . 

1  o 

Differentiation  of  this  expression  gives 


(5.7) 


x„  =  0  ,  <a  <  X,  <  X  : 
2  1  P 


)X  F' (x  /n) 

1  o 


(5.8) 


where  x,  =  x  defines  the  elastic-plastic  boundary. 

Ip 

At  the  elastic  plastic  boundary  the  deformation  gradients  are  and  , 


i.e.  , 


(5.9) 


It  can  be  shown  that  both  '!X,/;‘X,  and  3x,/3X„  are  continuous  across  x,  =  x  .  The 

112  2  Ip 

uniformity  of  and  in  the  intermediate  zone  also  implies  that  F^  and 

F^  depend  on  k  and  the  elastic  constants  only,  i.e.  F^  and  F^  are  independent 
of  a.  Equations  (5.8)  and  (5.9)  then  give  the  relation 


!  X  -  v''3a/< 

at  X  =  X  :  F-(  - )=  ^  , 

1  P  \  n  /  F,  ’ 


(5.10) 


where  (5.4)  has  also  been  used.  The  argument  of  F'(-)  depends  on  a,  but  the 

right-hand  side  of  the  equation  is  independent  of  a.  It  then  follows  that 

F(-)  is  either  linear  in  its  argument,  or  the  argument  is  constant.  The 

latter  case  would  be  x  -  ‘3a/-  =  constant.  Since  the  plastic  zone  should 

P 

vanish  as  a  1  0,  the  constant  must  be  zero,  implying  x^  <  a ,  which  is 
unacceptable.  Hence  (5.10)  can  only  give  a  linear  relation  for  F(-),  The 


integration  constant  has  been  set  equal  to  zero,  because  x  -  0  as  a  .  0. 

P 

If  we  would  take  a  as  the  argument  we  can  write 


5 .  Matching  of  Solutions 

It  may  be  assumed  that  the  particle  velocities  as  well  as  their 
derivatives  with  respect  to  are  continuous  at  =  <a  on  the  crack  line. 
It  then  follows  from  (3.18a)  and  (4.3b)  that  for  Ka  <  x.  <  x  ,  x„  =  0 

1  p  2 


V  =  —  a 

1  K 


(5.1) 


Equations  (4.4b)  and  (3.18c)  yield 


3  2 

K  a 


(5.2) 


The  function  C(t)  follows  by  matching  (4.6)  to  (3.18b).  Equation  (4.6) 
may  then  be  written  as 


_  /3a  1+6  2/3  a  6 

""2,2  “  3  2  1-6  ^1  ■  2  l-'S 

K  a  K  a 


(5.3) 


With  (5.1)  and  the  assumption  that  a  =  0  at  t  =  0,  we  may  rewrite 


(4,9)  as 


X  =  x_  - 
o  1 


(5.4) 


and  hence  (4.8)  becomes 


X,  =  X  (Xt  -  v'3a/<) 

1  o  i 

Continuity  of  X^  at  x^  =  <a  then  yields  the  identity 


(5.5) 


F(a)  =  X  (na)  ,  n  =  k  -  /3/<  , 
o 


(5 .6a,b) 


where  (3.28)  and  (5.5)  have  been  used.  Equations  (5.6a)  and  (5.5)  now 


imply 
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Here  F  and  G  are  arbitrary  functions,'  and  v„  »(s,x  )  follows  from  (4.6) 

2 , 2  o 

and  (4.9)  as 


(4.13) 


From  (4.11)  we  conclude  that  on  the  crack  line 


h(t,x  ) 


-r^  =  F(x  )  e 
3X2  o 


(4.14) 


In  the  next  section  the  functions  v, (t),  v,  -„(t),  x  (x, ,t),  C(t),  F(x  ) ,  X  (x  ) 

1  1,22  o  1  000 

and  G(x^)  will  be  determined  by  the  use  of  appropriate  matching  conditions. 


One  of  these  functions  is  actually  redundant,  but  the  form  (4.11)  is 
convenient  for  later  use. 
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'^2,2  1-6  '^1,22^’^^’^!  ’ 


(4.6) 


where  6  is  defined  by  (3.11),  and  C(t)  is  an  arbitrary  function. 

With  these  results,  one  can  relate  X  to  x.  Indeed,  since  X  h  0,  we 


3X  3X 

IT  ^  ^  ^  °  * 


(4.7) 


where  we  have  used  that  =  0  on  the  crack  line.  Equation  (4.7)  yields 


X.  =  X^(x^)  , 
1  0  0 


(4.8) 


where  X^  is  an  arbitrary  function  and  is  defined  by 


I. 

X  =  X,  -  I  V, 

o  1  1 


(s)ds . 


(4.9) 


To  determine  an  expression  for  Sx^/oX^  on  the  crack  line,  we  use 
(3/3x2)(X2)  =  0  ,  to  obtain 


Tt\3x„/  '^l  3x  \3x„/  '^2,2  9x„  ’ 

2  1^  2 


(4.10) 


where  we  have  used  that  v^  =  0  and  X^  H  0  on  the  crack  line.  The  solution 


to  (4.10)  is 


jX  -h(t,x  ) 

_ 2  _  I  o 

jx  F(x  )  ® 

2  o 


where 


(4.11) 


-V  ■  / 


G(x  ) 
o 


v_  ^(s,x  )ds. 
2,2  o 


(4.12) 
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4.  Crack-Line  Deformation  in  the  Zone  of  Uniform  Stress 

Next  we  investigate  the  crack-line  displacement  gradients  away  from  the 
tip,  in  the  region  where  the  stresses  are  defined  by  (2.10),  i.e., 

=  k  and  =  2k.  We  essentially  follow  the  procedure  proposed  by 
Achenbach  and  Li  (1984) . 

Substitution  of  (2.10)  into  (1.8)  and  (1.3)  gives  near  the  crack-line: 


/  0  -k  0 

i-5  <''1.2  (-•'  0  “ 

^0  0  0 


(4.1a) 


/  0  0  0 

T'  =  0  k  0 

\  0  0  -k 


(4.1b) 


Equation  (1.13)  then  yields 


-k 

• 

1  ° 

0 

°\ 

-k 

0 

0 

+  A 

0 

k 

0 

\0 

0 

0/ 

Vo 

0 

-k  / 

(4.2) 


Since  H  0,  we  have  on  the  crack  line 


v^  ^  =  0  ,  or  v^  =  Vj(t) 


(4 . 3a , b) 


Equation  (4.3a)  also  holds  near  the  crack  line,  and  hence  we  may  write  for 
=  0,  <  a  < 


'’1,221  "  '’1,22  ■  ''l,22<«- 

In  the  same  manner  (3/jx2)Dj^,,  yields  for  x^  =  0,  <  a  <  x^; 


(4.4a,b) 


(1+d)  '^2,21  ’ 


(4.5) 
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Equation  (3.25)  implies 


R  -  F(a.‘8<”>)  , 


where  the  functional  form  F(*)  will  be  determined  in  the  sequel.  Since 
the  displacement  is  radial  it  follows  from  (3.27)  that 


X  =^F(ae-2(P^ 
a  r 


where  g(p)  is  defined  by  (3.26).  Differentiation  of  (3.28)  yields 


3X 

I 

3x 


-  ■  7  (‘as  -  F'(ar8(‘>>)  f 

g  r  \a6  ^2  /  ^2  P- 


-g(p) 


f(p)  ■ 


On  the  crack  line  at  r  =  Ka  we  have 


3x 


1  K-f (k) 


3x, 


Ka 


3x, 


3X^  F'(a)  ’  dX^  F(a)  ’  3X2 


3x 

3X 


^  =  0 


(3.27) 


(3.28) 


(3.29) 


(3.30a,b,c) 
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Equation  v,3-13)  or  (3.15),  together  with  (2.7a)  and  (3 . 4a)  .determine  the 
velocity  field  near  the  tip. 

We  next  determine  the  relation  between  x  and  X.  To  this  end,  we  note 


3R  ,  3R  - 

—  +  V  -  =  0 

jt  a  3x 


(3.19) 


holds,  where 


R  =  (XX  )^. 
a  a 


(3.20) 


Equation  (3.19)  implies  that  R  =  const  on 


dt  ^  '^a  ’ 


(3.21) 


(3.22) 


where  we  have  used  (3.4a,b)  and  changed  the  independent  variable  from  t  to 
a.  It  is  more  convenient  to  use  p  =  r/a  instead  of  a  in  (3.22).  Noting  the 
relation 


da  _  1_  ^  _  r_ 

dp  p  dp  ; 


(3.23) 


we  transform  (3.22)  into 


f(p)  r 
p[p  -  f(p)] 


(3.24) 


which  integrates  to 


r/r,(R)  =  P  e  . 

where  r^(R)  is  an  arbitrary  function  of  R,  and 


(3.25) 


dC 

c-f(c)  ’ 


(3.2h) 
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where  (3.4a)  has  been  used.  The  exact  solution  to  (3.10)  and  (3.12)  can 
easily  be  written  out: 


f(a)) 


-rr/2/3  oj/3/2 

e _  e _ 

(2//3)^  (sinw)'"^ 


•l-sln(<.-V3)]'^/l  _ 
cos  (u)-it/3)  j  I 


where 


•n-/3 


-  / 


-5/3 

e  d5 


cos^  (^--n /3)  [  l-sin(5-n/3)  ]  *^ 


l.'hen  6  is  small,  which  is  usually  the  case,  this  solution  reduces  to 


f(w)  = 


-^/2/3  w  ►'3/2 
e _  e _ 

(2/k'3)^  (sino))'"^ 


(3.13) 


(3.14) 


(3.15) 


This  expression  satisfies  (3.12).  At  the  boundary  of  the  near-tip  zone, 
i.e.  for  oj  =  Tr/6  and  r  =  <a  we  find 


f (co  =  Tr/6)  H  f(p=  k)  =  /3/k  ,  (3.16) 

where  k  is  defined  by  (2.9). 

Equation  (3.15)  together  with  (2.7a)  and  (3.4a)  determines  the  velocity 
field  in  the  near-tip  zone.  For  future  reference  we  list  the  following 
results,  which  follow  from  (3.4a,b)  : 


,  =  A  f  (p) 

2,2  a  p 


'1,22  , 

a 


^  ^  f '  (p)?-f  C'") 
„  2 


Using  the  result  (3.16)  we  find  at  r  =  <a  : 


V  =  -  /3a 

1  <  ’  2,2  I 

K  a 


1,22 


*'  3a 


<'a^ 


(3.17a,b) 


(3.18a,b,c) 


Since  W  vanishes  identically,  we  have 


T  =  T  =  a 


h  T 


+  V 


a 

3X 


a 

By  the  use  of  (3.5),  (3.4  c-f)  and  (1.13)  we  obtain  equations  for  3v 
and  V  /r.  Subsequent  use  of  (3.4)  reduces  these  equations  to 


/3r 

r 


[f(0)-p!(-T;  +  2TJ)  +  x(|  . 

In  (3.6)  -  (3.7),  a  prime  denotes  a  derivative  with  respect  to  p. 
Elimination  of  A  from  (3.6)  and  (3.7)  yields  after  some  additional 
manipulation 


2T  -T„  .  ,  T'T  -T'T. 

cl/  \  r  0  f(p)  ,1  0  r  r  6  rc/»\  i  n 

f  (P)  -  2t  _T  p  2u  2T--T  tf(p)-p]  0. 
0  r  0  r 


According  to  (2.7a),  p  and  w  are  related  by 
tt/2/3  -uiv'3/2 

P  ^ 

P(^)  =  — 


(2//3)'''"  (sinoo)^" 


Elimination  of  p  from  (3.8)  by  the  use  of  (3.9)  then  yields 


duj 


cos  (aj+T/3) 
sinw 


cos  ((x)-tt/3) 


f  = 


<5p  (gj) 


cos(u)-7t/3) 


where 


6  =  k/2u  . 


The  boundary  condition  (3.3)  on  r  =  a  reduces  to 
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3 .  Crack-Line  Deformation  In  the  Near-Tip  Zone 
In  this  Section  the  displacement  gradients  are  determined  on  the  crack 
line  in  the  near-tip  zone  defined  by  a  <  <  Ka.  First  we  analyze  the 

velocity  vector  v  by  integrating  (1.13)  with  an  appropriate  condition  at 
r  =  a.  Since  we  have  assumed  that  near  the  crack  line  the  deformed  crack- 
tip  boundary  may  be  approximated  by  a  circular  segment  of  radius  a,  we  have 


XX  =  a^. 
a  a 


The  components  of  the  particle  velocity  then  satisfy 


(3.1) 


V  X  =  aa  ,  on  I x|  =  a  ^  (3.2) 

but  are  arbitrary  otherwise.  To  simplify  the  analysis  we  assume  that  v  is 
radial ,  1 .e. , 

V  =  —  X  on  |x|  =  a  .  (3.3) 

Then,  the  velocity  field  inside  the  material  also  becomes  radial  due  to  the 
hyperbolicity  of  the  governing  operation.  This  justifies  the  introduction  of 
the  following  forms: 

V  =  a  f(p)  ,  V  =  0  ,  (3.4a,b) 

r  u 

T^  =  T^(p)  ,  Tg  =  Tq(p)  ,  T^q=  0  ,  (3.4c,d.e) 

A  =  -  \(p)  , 

a 

where  p  =  r/a  ,  r  =  /(x  x  ) ,  and  f(p)  and  A (p )  are  functions  of  p,  and 

a  a 

T  (p)  and  T^(p)  are  the  stress  components  computed  from  (2.7b,c). 
r  y 


reason  we  have  selected  the  Thomason  field  for  the  intermediate  region. 

As  discussed  earlier,  however,  the  apex  angle  a  of  uniform  field  has  to  be  small. 
It  can  be  shown  that  the  point  A',  where  the  solution  (2.7)  is  connected  to 
a  uniform  field,  has  an  coordinate  smaller  than  Ka.  Since  the  system 
of  Fig.  4  reduces  to  that  of  Fig.  5  in  the  limit  of  a  4-  0  we  may  assume 
that  the  x^^  coordinate  of  A*  can  be  approximated  as  Ka  for  small  a. 

Accordingly,  one  may  say  that  the  stress  near  the  crack  line  beyond  A'  is 
uniform  and  approximately  equal  to 

=  k  1^2  =  2k  for  x^  >  Ka,  IX2I  ~  small.  (2.10) 

In  summary,  we  will  use  (2.7)  for  a  <  x^  <  Ka,  and  (2.10)  for 

Ka  <  X,  <  X  ,  where  x  is  the  x,  coordinate  of  the  elastic-plastic  boundary 
1  p  p  1 


on  the  crack  line. 


analytical  methods.  Consistent  with  the  approximate  nature  of  the  present 
analysis,  we  will  use  the  simple  'extension'  stated  above. 

To  discuss  the  matching  of  the  stress  fields  more  closely,  we  examine 
in  a  qualitative  manner  how  the  characteristics  of  the  fields  away  from  the 
tip  can  be  combined  with  those  of  the  near-tip  region.  Figures  4  and  5 
show  such  combined  systems  of  characteristics  for  the  Hutchinson  and  Thomason 
fields,  respectively.  To  the  left  of  the  characteristics  merging  at  points  A 
and  A'  the  stress  field  given  by  (2.7)  is  assumed  to  develop,  while  we 
assume  Hutchinson's  and  Thomason's  field  to  the  right  in  Figs.  4  and  5, 
respectively.  In  Fig.  4  the  crack  line  is  characteristic  to  the  right  of 
point  A.  Hence  the  upper  and  lower  characteristics  merge  at  A  without  slope 
discontinuity,  thus  forming  a  cusp  at  point  A  (which,  it  may  be  noted,  has 
an  x^  coordinate  of  <a).  From  this  observation,  we  may  infer  (and  can  prove) 
that  the  stress  at  point  A  has  essentially  the  same  mathematical  structure  as 
the  Hutchinson  field  (Fig.  2)  has  at  its  origin  (the  crack  tip).  This  con¬ 
sideration  and  the  analysis  by  Achenbach  and  Li  (1984)  then  suggest  that  the 
system  of  characteristic  lines  shown  in  Fig.  4  lead  to  an  anomalous 
deformation  having  a  singularity  at  A.  The  details  will  not  be  given  here, 
but  they  can  be  shown  by  using  the  same  method  as  reported  in  the  sequel. 
Since  a  deformation  singularity  at  A  would  be  unacceptable  we  are  left  with 
the  system  of  Fig.  5.  Due  to  the  presence  of  the  uniform  field  to  the  left 
of  the  point  A'  the  two  characteristics  merging  at  A'  will  not  form  a  cusp, 
and,  hence,  will  not  cause  an  unacceptable  singularity  there.  For  that 
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shown  in  Fig.  3.  The  essential  difference  with  Fig.  2  is  in  the  presence 
of  a  uniform  field  in  a  wedge-like  region  occupying  0  <  6  <  a.  Away  from 
this  wedge-like  region  the  two  fields  are  essentially  the  same.  Indeed, 
the  field  of  Fig.  3  reduces  to  that  of  Hutchinson  as  a  tends  to  zero.  For 
reasons  which  will  be  discussed  shortly,  we  will  adopt  this  solution  for  the 
intermediate  region,  but  we  will  assume  a  small  value  for  a.  Such  a  choice  of 
a  will  make  the  overall  stress  field  close  to  Hutchinson's  prediction  which, 
as  noted  earlier  appears  to  agree  with  other  calculations,  Hutchinson  (1968). 

In  fact,  the  uniform  field  in  the  wedge-like  region  differs  only  slightly 
from  =  k,  ~  2k  when  a  is  small. 

2 . 3  Matching  of  stress  fields 

In  order  to  obtain  a  complete  picture  of  the  crack-line  stresses,  we  have 
to  match  the  stress  fields  in  the  near-tip  zone  to  the  ones  in  the  intermediate 
region.  The  simplest  way  of  achieving  this  is  to  extend  the  near  tip  solution 
(2.7)  until  the  condition  T^^  =  k,  =  2k  is  met.  Equation  (2.7a)  shows 
that  this  condition  is  satisfied  at 

-  =  K  =  (/3  s  2.07  (u)  =  tt/6)  .  (2.9) 

One  would,  of  course,  expect  that  T  will  develop  dependence  on  9  in  the  transition 
from  the  essentially  axisymmetric  field  very  close  to  the  tip  to  the  non- 
axisymmetric  one  away  from  the  tip.  It  is  beyond  the  scope  of  the  present 
paper  to  analyze  this  transition  rigorously,  because  this  would  require 
more  information  on  the  near  tip  deformation  than  can  be  obtained  by 
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Figure  6  shows  the  elastic-plastic  boundary  in  the  x^-a  plane.  The  "a" 

coordinate  of  the  intersection  of  this  curve  and  the  line  defined  by  (5.4)  on 

which  X  =  constant,  defines  a  function  of  x  •  We  shall  choose  this  function 
o  o 

*  it 

for  G  (x^)  i^  (5.14).  Since  this  construction  of  G  yields 

G*(x  )  =  a  (5.18) 

o 

on  the  elastic-plastic  boundary,  we  find  from  (5.9)  and  (5.14)  that 

F(x^)  =  F^  .  (5.19) 

Continuity  of  at  x^  =  xa  gives  by  the  use  of  (3.30b),  (5.11) 


and  (5.14) 


nP~l^  e^n  h  ^  /  na  ^qna/G  (na) 

^2  VG*(na)/ 


(5.20) 


It  follows  from  (5.20)  that 


G*(na) 

ha 


=  y 


(5.21) 


where  y  is  the  root  of  the  equation 


pP  =  Y~P 

^2 


Hence,  we  have 


G  (x  )  =  Y  X  , 
o  o 


Equations  (5.14),  (5.19)  and  (5.23)  then  determine  in  ca 


(5.22) 


(5.23) 


clX„  .  p 


-qx  ( - ) 

o  a  -fx 
e  o 


X,  <  X  as 

1  p 


(5.24) 


6.  Results 


On  the  crack  line  the  near-tip  region  corresponds  to  a  ^  ^  Ka, 

where  k  is  given  by  (2.9):  <  =  2.07.  Equation  (3.28)  relates  to  by 

X,  .  FCae-Sf-*)  .  ^  .  (6.1) 

^  ^1 

where  (5.11)  has  been  used  and  g(p)  is  defined  by  (3.26).  Here  p  should  be 
interpreted  as  p  =  x^/a,  and  n  is  given  by  (5.6b).  For  x^  =  Ka,  (6.1)  yields 
=  pa/F^.  At  x^  =  a  we  have  p  =  1  and  (6.1)  yields  X^  =  0,  because 
f(l)  =  1,  and  hence  g(p)  -  [f'(l)-l]  ^iln(p-l)  as  p  'I-  1,  which  implies 
exp[-g(p)]  -  (p-1)^^^  as  p  4-  1,  because  (3.9)  and  (3.10)  give  f'(l)  =  -  h- 
The  crack-line  stresses  in  the  near-tip  zone  are  given  by  (3.4c,d)  and 
(2.7b,c).  Thus, 


0  ^  Xj^  _<  na/F^,  or  a  _<  xa: 

T^^  =  2kcos(uJ  +  ti/6)  ,  T22  2kcos  (w-tt/O)  , 

where  't/6  w  _<  :r/3  ,  and  according  to  (2.7a)  x^  and  co  are  related  by 

/l^^  /3  v'3(t/3-(d) 

2  sini.) 


The  deformation  gradients  on  the  crack  line  follow  from  (3.29): 


iX., 


(6 . 2a ,b) 


(6.3) 


V 


(6.4) 


At  =  a,  Sx^/SX^  is  singular,  but  3Xj^/3X^  vanishes. 

The  intermediate  region  is  defined  by  o  ^  ^  x^.  The  corresponding 

range  of  X,  follows  from  (5.12)  as  na/F,  <  X^  <  (x  -/3a/ic)/F,  .  The  crack- 
i  1  —  1  —  p  1 

line  stresses  in  the  intermediate  region  are  given  by  (2.10).  Thus 


pa 


/F,  <  X^  <  (x  -t/3/K)/F,  ,  or  Ka  <  X,  <  X  : 

1—1  —  p  1  —  1—  p 


Til  =  k  .  ^22  = 


From  (5.12)  we  conclude 


while  Sx^/SX^  is  given  by  (5.24). 


(6.6a,b) 


(6.7) 


The  computation  of  numerical  results  requires  the  evaluation  of  the 
integral  in  (3.26),  where  f(p)  follows  from  either  (3.13)  or  (3.15)  and 
(2.7a).  The  precise  evaluation  of  the  integral,  which  is  not  difficult 
numerically,  may  not  be  worth  the  effort  considering  the  approximate  nature 
of  the  present  analysis.  Hence,  we  will  simplify  the  analysis  by  using 
the  approximate  formula  (3.15)  for  f.  In  addition  it  is  found  that  the 

A 

Hermitian  interpolation  f  of  f 


=  -L(p-l)(p^+C^p 


+  C^) 


+  P 


(6.8) 


is  a  sufficiently  accurate  approximation  for  f,  where 
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=  - 


2{k^  -  4k^  -  (l-3/3)<  +/3} 
-  5k  +  4/3 


.  3k **  -  7k ^  +  6/3k  -  2/3 
K^  -  5k  +  4/3 

3 

2(1+C^+C2)  ■ 


(6.9a) 


(6.9b) 


(6.9c) 


It  is  noted  that  this  approximation  will  predict  the  correct  singularity 

for  deformation  gradients  at  the  crack  tip  because  of  the  use  of  Hermitian 

interpolation.  Indeed,  as  was  seen  earlier,  the  singularity  of  9x/3X  is 

determined  by  the  derivative  of  f  at  p  =  1.  Also,  it  is  evident,  from  (3.8), 

that  the  singularity  of  3x/3X  does  not  depend  on  6.  This  observation  serves 

* 

as  another  justification  for  the  use  of  (3.15).  With  f  replacing  f  in 
(  3. 2»-/) ,  we  obtain 


3  2 


,  2/3 

s 

/  K-W  \  O 

(  °  ) 

(p-l ) 

\  P-UJ  / 

\  / 

0 

(6.10) 


(6.11) 


where 


-C^+i  (4C2-cp' 


6  =  -  -  +  i 

o  3 


^  3  (4C2-C^f^ 


(6.12a,b) 


We  also  have 


/  1  v2/3 

3  1 

/  P -w  \  0 

(  °  ) 

^<-1  ) 

\  <“/A)  / 

O  ! 

(6.13) 
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Hence  we  have 


Further  simplifications  may  be  introduced  when  6  =  k,/2p  <<  1.  Then 


(6.14) 


F  „  =  1  +  0(6)  £  1  ,  p  s  3/^“  s  0.16,  q  =  /3/<^  s  0.19. 


(6.15a,b,c) 


With  these  constants,  (5.22)  gives 


Y  =  0.321. 


(6.16) 


This  result,  (5.4),  (5.23)  and  (5.18)  yield  x  =  3.95a  and 

P 

from  (5.12),  we  see  that  the  material  in  the  range  0  <  <  3.12a  will  be 

in  plastic  zone.  Plots  of  the  deformation  gradients  on  the  crack  line 
for  6=0  are  given  in  Fig.  7. 

Finally,  we  can  make  an  estimate  of  a  in  terms  of  the  stress  intensity 
factor  Kj  of  linear  elastic  fracture  mechanics  by  equating  the  plastic  zone 
sizes  predicted  by  the  present  theory  and  existing  small  strain  theories. 
Small  strain  theories  predict  that  the  length  of  the  plastic  zone  x^  is 


given  by 


X  =  c  yI/V}  , 
P  I 


(6.17) 


where  c  =  tt/24  =  0.13  according  to  the  Dugdale  model,  and  c  =  2/2/9it  s  0.10, 
according  to  the  analysis  Achenbach  and  Dunayevsky  (1984).  By  equating  (6.17) 


to  X  -a  =  2.95a  and  X  =  3.12a  respectivelv ,  we  obtain 
P  P 

~  c  I  J  _  C  I 

S  -  <5  Q  c  3nQ  s  -  • 

2.95  J.12  ^2 


(6.18a,b) 
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Fig.  2  Characteristic  curves  for  Hutchinson's  field.  Point 
0  corresponds  to  the  crack,  tip. 


Fig.  3  Characteristic  curves  for  Thomason's  field.  Point  0 
corresponds  to  the  crack  tip. 
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Fig.  4  Characteristic  curves  for  crack-tip  zone  and  intermediate 
region,  with  Hutchinson's  field  in  intermediate  region. 
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Fig.  5  Characteristic  curves  for  crack-tip  zone  and  intermediate 
region,  with  Thomason's  field  in  intermediate  region 
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